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Abstract: The transport of the nuclear cosmic ray flux in the atmosphere is studied and the atmospheric 
corrections to be applied to the measurements are calculated. The contribution of the calculated correc- 
tions to the accuracy of the experimental results are discussed and evaluated over the kinetic energy range 
10-10 3 GeV/n. The Boron (B) and Carbon (C) elements system is used as a test case. It is shown that 
the required corrections become largely dominant at the highest energies investigated. The results are 
discussed. 



Introduction 

The direct measurement of nuclear cosmic 
rays (CR) by balloon borne experiments like 
CREAM [9] is a precious tool to understand the 
cosmic rays propagation process in the galaxy. 
Among the different observables, the secondary- 
to-primary ratios of nuclear CR like B/C play a pre- 
dominant role because they measure directly the 
thickness of matter crossed by CR and therefore 
tell us about their confinement in the galaxy. In ad- 
dition, the CR flux entering the atmosphere under- 
goes nuclear fragmentation by interaction with the 
atmospheric nuclei, producing more secondaries 
and therefore enhancing the measured secondary- 
to-primary ratios. 

At high energy, the confinement of cosmic rays in 
the galaxy becomes less and less efficient. The 
thickness of matter crossed in the galaxy is there- 
fore decreasing with energy whereas the thickness 
of matter crossed in the atmosphere overhead the 
detector remains constant with energy. The pro- 
duction of secondaries in the atmosphere is then 
becoming dominant at high energy. 

The purpose of this work is to propose a framework 
to compute the evolution of the nuclear cosmic ray 
composition in the atmosphere, taking into account 



the absorption and fragmentation processes. It can 
be used to reconstruct the top of atmosphere (TOA) 
flux from the value measured by a detector in the 
atmosphere. The effects of statistical and system- 
atic errors on this reconstruction are also carefully 
studied. 

Transport model 

The weighted slab model (WSM) provides a partic- 
ularly appropriate framework for the calculations 
since in this approach, the variable of the transport 
equation is the amount of matter crossed by the 
particles (grammage), see Ref. [4] and references 
of this paper for details on the WSM. In this model 
the transport equation can be written at a given ki- 
netic energy per nucleon E k , as: 

dN z (x,E k ) = y^ gr(g) x ( 
dx Z-j mT 

I - <r iT (E k )Ni(x,E k ) + J2 aJ iT N j( x i E k) 

I 3>i 

where Ni is the abundance of nuclear element i at 
atmospheric depth x, tut being the target mass, 
am the total reaction cross section and aj the 
fragmentation cross section from a nucleus j to 
a (lighter, i < j) nucleus i. The TOA flux from 
ref [1 1] were solar modulated and used for the ini- 
tial conditions A^(0). The sum over T corresponds 
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to the various components of the atmosphere, with 
mass rriT and fraction R T (x) = pT ^f\ which 

1 1 V J ptot{x) 

does not change significantly for the altitudes be- 
tween ground and 200 km for the three main con- 
stituents [3]. The numerical approach of the prob- 
lem has to be handled with care. The inversion of 
equation 1 where the TOA flux has to be computed 
from the measured values leads to numerical dif- 
ficulties and lengthy calculations if the direct inte- 
gration method is used. In this work, the numeri- 
cal calculations were instead, performed using the 
simpler, easier to handle, matrix formulation of the 
problem, in which the inversion is easy to achieve. 

In this framework, the transport equation 1 can be 
expressed as: 

AN{X Ax Ek) =S(x,E k )N{x,E k ) (2) 

where N(x, E k ) is the vector containing all the el- 
emental abundances Ni(x,E k ) of the considered 
CR flux for a traversed grammage of x g/cm 2 in 
the atmosphere. S is the transformation matrix. It 
is a triangular matrix, with the diagonal elements 
of S corresponding to the nuclear absorption and 
the other elements to the production of secondary 
nuclei. The solution of 2 is given by: 

N(x,E k )=R(x,E k )N(0,E k ) (3) 

where R(x,E k ) = exp (S(x,E k )) is the transfer 
(transport) matrix. 

To compute the S and R matrix, the parametriza- 
tion from [8] was used for the total reaction cross 
sections a it, while the fit formula from [10] was 
used for the fragmentation cross sections o\ T . 
In this latter case however, the formula applied 
mainly to H and He targets. Its application was ex- 
tended here to larger masses using a factorization 
approximation [6]. 

The atmosphere model from [3] was used in these 
calculations. 

Results 

Figure 1 shows the non-propagated and propagated 
cosmic ray fluxes from a calculation for a spectral 
index of the diffusion coefficient 6 = 0.3, versus 
the energy per nucleon and the correspondant B/C 
ratio for different values of the spectral index of 



the diffusion coefficient 5 = 0.3, 0.46, 0.6, 0.7, and 
0.85. It can be seen that at the highest energy con- 
sidered, the flux values at TOA (solid lines) and 
at the balloon detection altitude after crossing 5 
g/cm 2 of atmosphere (dashed lines), differ by ap- 
proximately a factor of two for small values of 5, 
and by more than one order of magnitude for the 
largest values. This is easy to understand since for 
larger S the galactic contribution at high energy be- 
comes dominated by the atmospheric contribution. 
The asymptotic limit of this behavior can be ob- 
served with the propagated ratio tending to flatten 
out at high energy above 1 TeV/n for large values 
of S where the galactic secondaries production be- 
comes negligible. 

This illustrative example shows how critical are the 
corrections to be applied to the measured flux val- 
ues in balloon experiments, and thus how impor- 
tant is a careful study of CR transport in the atmo- 
sphere for a reliable results in TOA flux evaluations 
from atmospheric measurements. In the high en- 
ergy region, the measured raw B/C value appears 
to be about one order of magnitude larger than the 
TOA flux value to be extracted from it. 

The flux at TOA can be reconstructed from the flux 
measured at a given thickness x by inverting the 
equation 3: 

N(0,E k ) = R(x,E k )- 1 N(x,E k ) 

But in the inversion process, distortions may be 
generated by the random (statistical) fluctuations 
of the experimentally measured fluxes N(x,E k ). 
Correcting for these distortions is equivalent to 
the unfolding [2] of a measured spectrum from 
the detector response, the role of the latter being 
played here by the transfer function (inverted ma- 
trix). This effect has been shown to be negligi- 
ble for the grammage ( 5 g/cm 2 ) considered in this 
study [7]. 

The uncertainties induced on the flux calculations 
by the experimental uncertainties on the nuclear 
cross-sections have been estimated by the follow- 
ing way for three values of the spectral index of the 
diffusion coefficient 5=03, 0.6 and 0.85, covering 
the range of realistic values [5, 1]: 

A sample of one hundred transfer matrix 
R err (x,E k ) was generated by adding ran- 
domly a ±5% systematic error to the total reaction 
cross sections (diagonal elements), and ±10% 
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B/C were searched in the 100 values calculated 
with the error-added matrices, and used as an es- 
timate of the upper and lower limits of the uncer- 
tainties induced by the systematic errors on cross 
sections. Figure 2 shows the results from these cal- 
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Figure 2: The filled areas show the systematic er- 
ror region obtained by the method described in the 
text, for the three different values of S = 0.3, 0.6, 
0.85. The statistical error appear as the error bars 
on the figure. They were calculated for 100 days 
of effective counting and for 15 energy bins per 
decade of energy. 



Figure 1 : Results of the transport calculations for 
C and B nuclei. Top: Differential fluxes with 5 = 
0.3 at TOA (solid lines) and for a crossed matter 
thickness of 5 g/cm 2 (dashed lines) for 12 C (upper 
curves) and n B (lower curves). Bottom: B/C ratio 
at TOA (solid lines) and after propagation (dashed 
lines) for 6: 0.3, 0.46, 0.6, 0.7, 0.85, from top to 
bottom respectively. All curves as a function of the 
kinetic energy per nucleon. 

to the fragmentation cross sections (off diagonal 
elements) in R(x, E k ). The measured fluxes were 
calculated using the error added matrix 

N^(x,E k ) = R e ™(x,E k )N(0,E k ), 

and the TOA fluxes were reconstructed by the in- 
version procedure using the nominal matrix R(x): 

7V err (0,^) = R-\x,E k )N^(x,E k ), 

Then the B/C ratios were calculated for each en- 
ergy Ek, and the minimal and maximal values of 



culations for three values of the spectral index of 
the diffusion coefficient 5=0.3, 0.6 and 0.85, cov- 
ering the range of realistic values [5, 1]. The val- 
ues obtained from the inversion of R(x) are taken 
as central values. These results show that as ex- 
pected from fig 1, for 6 = 0.3 the systematic er- 
ror is increasing with energy, but remains relatively 
small over the whole energy range, in contrast with 
the S = 0.85 result, where it becomes large, ex- 
tending over almost a factor 2 at high energy. As 
already mentioned, this can be easily understood 
since for S — 0.85 more secondaries are produced 
in the Earth's atmosphere compared to the sec- 
ondary galactic flux compared with the S = 0.3 
productions (see figure 1). 

Another question to be addressed is what count- 
ing statistics is needed for the measurements to be 
dominated by systematics errors. A simple evalu- 
ation was made by assuming a detector acceptance 
of 1 m 2 sr, and taking 30 bins over the energy 
range 10-1000 GeV/n, and 100 days of measure- 
ment time. The obtained values are shown as error 
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bars on the B/C ratios on Fig. 2 for the three values 
of spectral indices S. It can be seen on the figure 
that the statistical error is dominant for high ener- 
gies (Ek > 200 GeV/n) for the assumed experi- 
mental condition. A more quantitative study com- 
parable to that reported in [1] remains to be per- 
formed. In this reference, a x 2 minimization pro- 
cedure on the simulated B/C data was performed, 
to evaluate the accuracy in the determination of S 
due to statistical errors, with all the propagation pa- 
rameters left free. The estimated statistical uncer- 
tainty on the determination of 5 was 10-15% for 
experiments like CREAM. 

Summary and conclusion 

The raw CR flux of nuclear elements measured in 
balloon experiments must undergo significant cor- 
rections for the effects of propagation in the atmo- 
sphere, to provide a reliable measurement of the 
TOA fluxes. These corrections become larger with 
the increasing energy and with increasing diffu- 
sion spectral index. Due to the uncertainties on the 
absorption and fragmentation cross sections, they 
dramatically affect the accuracy of the experimen- 
tal values at high energy, but the measurements 
keep a good part of their capacity to discriminate 
between the values of S predicted by the various 
models. 

Since the same fragmentation process takes place 
in the detectors, similar results can be expected, 
and the measured raw flux would have to be cor- 
rected by a similar procedure as reported in this 
paper to account for the flux propagation inside the 
detector material. The correction will of course de- 
pend on the architecture of the apparatus of the ex- 
periment and of the amount of matter upstream of 
the charge measurement detector(s). 
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